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Abstract 

The metastable states of a glass are counted by adding a weak pinning field 
which explicitly breaks the ergodicity. Their entropy, that is the logarithm 
of their number, is extensive in a range of temperatures Tq < T < Tc only, 
where Tq and Tq correspond to the ideal calorimetric and kinetic glass tran- 
sition temperatures respectively. An explicit self-consistent computation of 

the metastable states entropy for a non disordered model is given. 
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During the last decade, a great deal of work has been devoted to the understanding 
of the glass transition. Basically, two different approaches have been employed to tackle 
this problem. On the one hand, the density functional theory (DFT) Jl] is an equilibrium 
approach where the static density field obtained through the minimization of a free-energy 
functional becomes inhomogeneous below the structural glass transition temperature Tq. 
On the other hand, the mode coupling theory (MCT) is an off-equilibrium theory which 
self-consistently treats the microscopic dynamical correlation and response functions in the 
liquid phase. The onset of the glassy state is then due to a kinetic transition at a temperature 
Tc higher than Tq- The breaking of ergodicity at Tc coincides with the appearance of non 
vanishing density fluctuations in the long time limit, the system becoming partially frozen 
in metastables states with very large relaxation times. 

The connection between the dynamic and the static approaches can been made more 
transparent and rigorous in the special case of mean-field spin glasses |§f|]. It has indeed 
been found that the models exhibiting a discontinuous replica symmetry breaking || at the 
static transition temperature Tg present a dynamic transition at a higher temperature Tp 
below which equilibrium is never reached In addition, the temperature T D whose 

significance is a priori purely dynamical may be computed from the Gibbs partition function 
using the so-called "marginality condition" ||. This is related to the intuitive feeling that 
the dynamical transition should arise when the free-energy landscape becomes "rough" and 
includes many states separated by extensive barriers while at the same time the replica 
symmetry remains unbroken (To > T5) . Recently, it has been shown that these results may 
be extended to mean-field spin models with complicated but non-random interactions 
A glassy system, once trapped in a metastable state, is indeed partially frozen and its slow 
degrees of freedom may act as self-induced "quenched" constraints with respect to the fast 
ones m§- 

In this letter, we wish to investigate the metastable states in the equilibrium free-energy 
landscape in any glassy system with or without quenched disorder. More precisely, we shall 
define the entropy Sh s as the logarithm of the number of these "hidden" metastable states 
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and present a general scheme to compute this quantity. It will be shown that for mean-field 
disordered models Sh s is extensive in the range of temperatures T s < T < T D only, and 
gives also some important information about the structure of the states below T$. We shall 
argue that the same holds in the absence of disorder and that these two temperatures then 
correspond to the ideal calorimetric and kinetic temperatures Tq and Tq respectively. The 
replica formalism will then be adapted to the case of systems without disorder to derive 
their metastable states entropy Shs and the glassy correlation functions inside these states. 
An example of such a calculation will be given for a simple model. 

Let us start with a theory of a field <p(x) defined by a Hamiltonian H[<ft}. For simplicity, 
we shall use scalar notation though x spans a D-dimensional space and 4>(x) can be an 
iV-component field. The equilibrium Gibbs free-energy at the temperature T = jj is given 
by 

Fttf) = -hag J dfiix) e-WQ (1) 

For a usual ferromagnet, the emergence of a spontaneous magnetization at low temperature 
coincides with the breaking of the symmetry <p{x) — > —<j>{x) of the Hamiltonian. Below 
the Curie temperature, the physical decomposition of the Gibbs free-energy into two states 
of opposite magnetizations may be obtained by imposing a small (but finite) external field 
aligned along the up or down directions. In the case of disordered systems or glasses, there is 
however no a priori privileged direction towards which the field <p(x) will point once stuck in 
a metastable state. We can nevertheless choose a possible direction, given by another field 
<j(x), and compute the free-energy of our system when it is weakly pinned by this external 
quenched field 

[a, g, (3] = ~ log J d<j)(x) e-^M-f JM^)-H-)) 2 ( 2 ) 

where g > denotes the strength of the coupling. This free-energy (|2]) will be small when the 
external perturbing field cr(x) lies in a direction corresponding to the bottom of a well of the 
unperturbed free-energy ([!]). Therefore, we should be able to obtain useful information about 
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the free-energy landscape by scanning the entire space of the configurations a(x) to locate 
all the states in which the system can freeze after spontaneous ergodicity breaking (g — > 0). 
According to this intuitive idea, we now consider the field cx(x) as a thermalized variable 
with the "Hamiltonian" [a, g, (3\. The free-energy of the field a at inverse temperature 
(3m where m is a positive free parameter therefore reads 

FJm,(3) = lim - J- log / da(x) e - pmF ^' 9 ^ (3) 
g^o+ (3m J 

When the ratio m between the two temperatures is an integer, one can easily integrate ct(x) 
in eqn.(|3|) after having introduced m copies (b p {x) (p = l...m) of the original field to obtain 
the relation 

FJm,(3) = lim -—log / TT d<P P (x) e -^^M+IE P , A ^/^W^W (4) 
9^o+ (3m J 

where g pX = g{— — 5 pX ). Let us define two more quantities related to the field a : its 
internal energy W(m,(3) = and its entropy S(m,(3) = f3m 2 ^J^. Since the case 

m = 1 will be of particular interest, we shall use hereafter Fh s {/3) = W(m = 1,(3) and 
Sh s (P) = S(m = 1,(3) where hs stands for "hidden states". We stress that S(m,(3) and 
(3 2 ^-§p which are respectively the entropies of the fields a and are two distinct quantities 
with different physical meanings. 

When the pinning field <j(x) is thermalized at the same temperature as 4>(x), that is 
when m = 1, one sees from eqn.(^j) that = F a (m = 1, (3). The basic idea of this letter 

is to decompose F a into its energetic and entropic contributions to obtain 



S hs {(3) = (3 



F hs {(3) - FM 



(5) 



To get some insights on the significance of the above relation, we shall now turn to the par- 
ticular case of disordered mean-field systems. We shall see how it rigorously gives back some 
analytical results derived within the mean-field TAP and dynamical approaches [p|,|6||8||9|,|TTf . 
We shall then discuss the physical meaning of identity (||) for the general case of glassy 
systems. 



The free-energy F^fi) of a mean-field disordered system is a self-averaging quantity which 
may be computed using the replica trick H to end up with FJ(3) = lim Ext J r J{q ab }) . 

' n^O {q ab } 

q ab = y J dxcf) a (x)(f) b (x) are the overlaps between the n replicas (V = J dx is the volume 
of the system). Above the static transition temperature T$, the physical saddle-point of 
T$ is replica symmetric (RS) (for simplicity a spherical normalization q aa = 1 

has been assumed). If we now compute the free-energy @ of the field a by introducing n 
replicas and averaging over the quenched disorder, we obtain from eqn.@ the same free- 
energy functional T$ where the number of replicas (p ap (x) now equals n x m and with 
the additional term Sa=i 2JTa=i g px q ap,aX . This interaction term explicitly breaks the 
symmetry of permutations of the n x m replicas [jnj into n groups of m indistinguishable 
replicas. Consequently, even above T s , the simplest Ansatz one can resort to contains at 
least one step of replica symmetry breaking (RSB) which reads 

F a (m,f3) = Ext ^(g , ?i, m, /9) (6) 

90, 9i 

when g — > + . As a result of the introduction of the field a, we have obtained the usual 
one-step expression of the free-energy but without any optimization over the free parameter 
m which we can choose at our convenience (see \Y1\ for a related case where a acquires a 



simple geometrical interpretation). Let us now send m — > 1 while T > T$. The saddle-point 

equation over q becomes identical to the RS equation for q. Thus q = q and we find F a {m = 

dJ~, 

1,(3) = F^{0) = ExtXqJ 7 ^,^) as expected. Defining V(qi) = ( 3 ~g^( ( lo = 1^^ m = 
the optimization condition over the second overlap implies that q\ must be a stable local 
minimum of V. The entropy of a is then Sh s — V(<Zi). ^ turns out that V defined above 
is equal to the potential recently introduced in mean-field glasses P,|10| to compute the 
marginality temperature To at which the relaxational dynamics exhibits a drastic slow 
down as mentioned in the introduction. The typical behaviour of the entropy Sh S {P) 
is as follows. At high temperature, there is only the RS solution q\ = q, Shs = 0. At a 
given To, there appears a non trivial saddle-point qx (which shifts the free-energy (|5[) by an 
extensive amount of order gV ||15|| ) and the entropy Shs shows a first order jump. When T 



decreases, Shs goes down and vanishes at T$. Note that Sh s (T = T$) = is mathematically 
equivalent to the usual optimization condition of the one-step free-energy with respect to 
m. This is an important remark we shall comment on in the following. 

The marginality condition |||9|] tells us that T = Td is the temperature of the onset 
of an exponential number of such metastables states. Using the TAP formalism [JO]], one 
can indeed write the set of equations that the mean- field local magnetizations {M(x)} 
must satisfy. The paramagnetic solution M(x) = (for simplicity we suppose that q = 
(y J dxM(x)) 2 vanishes in the RS phase) gives the RS free-energy density as expected. 
Below Td the RS free-energy density is also achieved by the sum of an exponential number 
of solutions with a higher free-energy and an internal overlap qi — y J dxM(x) 2 ||. It has 
been recently verified in the particular case of the random orthogonal model || that the 
logarithm of the number of these solutions - which is sometimes called complexity or 



configurational entropy fl3,|14 - coincides with the entropy Sh s defined in this letter. 



The significance of formula (|5[) may now be discussed for glasses in general. At high 
temperature, the pinning due to a is not sufficient to break ergodicity : F hs (f3) = F^{(3) 
and Shs(P) = 0. Physically, there exists only one state. When the temperature goes down, 
there may appear some high barriers which separate an exponential number of metastable 
states of free-energies higher than the true F${f3). As long as the number of these states 
does not compensate for their small weights, they are not "seen" by the Gibbs partition 
function. At some temperature Td, their number e Shs becomes large enough to make up for 
the difference of free-energies and the identity (|5]) expresses this compensation mechanism. 
In a system with finite range interactions, the ergodicity breaking taking place at Td is not 
complete. The partial freezing of collective modes inside the metastable states makes the 
dynamics very slow but does not forbid some microscopic changes (the so-called activated 
processes) [[|. Thus, no discontinuties occur at T D for the thermodynamic quantities and 
it is reasonable to think that Td coincides with the ideal kinetic transition temperature Tc 
for glassy systems. Below Td, F^{f3) comes from the superposition of many states with high 
free-energies. The number e Shs of these hidden states decreases since their free-energy Fh s 



gets closer and closer to the true value F^. At a given temperature Tg, these states cease to 
be metastable since their free-energy F hs equals Fs, implying from (El) that the entropy Sh s 



is not extensive anymore |14J . The true thermodynamical transition therefore takes place at 
Ts which corresponds to the usual calorimetric glass transition temperature Tq. Below T$, 
formula (|5[) cannot hold any longer since it would predict a negative configurational entropy 
of metastable states having a lower free-energy than F^. Physically, one expects that freezing 
into a small (non exponential) number of states will still occur and that Sh s {T < Ts) = 0. 
As a consequence, the effective temperature of the states, that is of the field a, becomes 
higher than the temperature of the field <fr and m is determined through the condition that 
the hidden states' entropy vanishes. If one goes back to the case of mean-field disordered 
systems, one finds that m is such that the one-step free-energy derivative with respect to 
m equals zero. Therefore, we have obtained a simple physical interpretation of the usual 
optimization condition of T$ with respect to m in the RSB phase, whose meaning has 
always been far less clear than the optimization with respect to the overlaps q and q\ [Q. 
Furthermore, let us notice that if Shs is already non extensive at then Tp and Ts must 
coincide. This is what happens for systems with a continuous RSB transition || where the 
intermediate phase Ts < T < T D is skipped. Such systems seem therefore to exhibit a less 
generic behaviour |3||§]. 

Formula @ is a convenient starting point to compute Sh s in systems without quenched 
disorder. We begin with m uncoupled copies (p p (x) of the model. The matrix of the correla- 
tion functions of the global system is then a priori diagonal : G pX (x — y) = (<j) p (x)<j) x (y)) = 
S pX G(x — y). For simplicity we shall assume that in the liquid phase the average value of 
the field (4>(x)) is equal to zero. The onset of the glassy phase will be characterized by 
the appearance of metastable states in which the expectation value (<p(x)) does not vanish 
anymore and therefore by the emergence of non zero off-diagonal propagators 

G hs (x-y) = (<f ) (x)}( ( j ) (y)} = ( ( f ) p (x)<f> x (y)} (p ? A) (7) 

in the limit m — > 1 As discussed above, the diagonal correlation function G(x — y) is 



simply obtained from the calculation of F${f3) for a single system (m = 1). The computation 
of F a (m 1,(3) is more delicate since a perturbative expansion of the original Hamiltonian 
will never be able to generate an effective coupling between the replicas and to detect the 
first order glass transition. We shall now present a simple self-consistent calculation which 
permits us to obtain Sh s and Gh s (x — y). 

Let us consider a system of N particles at positions Xi (i = 1...JV) interacting through 
a two-body potential J2i<jV(\xi ~ x j\)- In the grand-canonical ensemble, this system is 
described by the Hamiltonian H[<f>] — hj dxdy — y)4>(y) — \i j dx exp i<f>(x) where 

H is the chemical potential and <p(x) the conjugated field to the density of particles at point 
x. To obtain an analytically tractable problem, we shall make two simplifications. First, 
while the previous general expression contains interactions at all orders in 0, we shall keep 
only the quartic vertex. Secondly, we shall consider a O(N) version of the above Hamiltonian 
and compute the free-energy F a (m,(3) using Bray's self-consistent screening approximation 
||17|| . The expansion parameter 1/N plays the role of the chemical potential \i. We do not 
expect this model to exhibit a glass transition for large N but the small N case may have 



an important physical interest, e.g. polymers which are related to N — > []T8| are known to 
have a glassy behaviour at low temperatures. Bray's calculation is exact to order 1/N and 
contains a partial resummation of an infinite class of diagrams |17| and is thus well defined 
in the whole range < iV < oo. With the Ansatz G pX (k) = 5 px G(k) + (1 - 5 px )G hs (k) 
||16|| , one obtains the free-energy F a (m,N) as an extremum over the set of all propagators 
G(k) and Gh s {k) jl9| . When m — > 1, the diagonal propagators are solutions of the implicit 
equations 

where H(k) = J dqG(k — q)G(q). Eqn.(||) is identical to Bray's result as expected [0]. The 
entropy of the metastable states then reads 



S hs (N) = Ext fdk 

G h3 (k)J 



s ( G hs (k) \ _1_J U hs (k) 



G(k) J N \1 + U(k) 



(9) 



where n^fc) = / dqGh s (k—q)Gh s {q) ands(x) = — x— log(l— x). The numerical resolution of 
the saddle-point equations steming from eqn.([8],|9]) in dimension D = 3 shows that a non zero 
set of propagators Ghs may appear when N becomes lower than a given Nc which depends on 
the bare propagator V{k). Despite different choices of V(k), we have always found Nc < 1. 
A more careful analysis of the equations would however be useful. In dimension D — 0, the 
equations for G and Ghs may be solved exactly and we find that there exists a first order 
transition at some small enough critical value of N. One can show that Nc is always lower 
than one. If for instance we choose the bare mass V = —0.2, Sh s is equal to zero when 
iV > Nc — 0.65 and jumps discontinuously to 0.31 at the transition with Gh s /G ~ 0.91. 
The entropy then decreases smoothly when N gets smaller and vanishes at Nc — 0.54 where 
Ghs/G ~ 0.96. We notice that the ratios Nc/Nc and Gh s /G are remarkably similar to the 
values of Tp/Ts and q\ which may be found in mean-field disordered models Though 
one must consider this result with caution due to the approximations made in its derivation, 
it seems that the free-energy landscape of the (<f) 2 ) 2 model may be complicated at small N 
(N < 1), even above the ferromagnetic transition temperature. This is strongly reminiscent 
of the random field Ising model which is recovered here if N = 1 and g keeps a finite value 

A deeper and more rigorous understanding of the mechanism of the ergodicity breaking 
occurring at Tc is still to be found. In this respect, one could try to transpose the TAP 
approach to non disordered models. Let us call T(M) = g J dxa(x)M(x) + F^[a,g,/3] the 
Legendre transform of the "Hamiltonian" of the field a. The free-energy F^,{J3) is equal to 
the minimum of the effective potential T at M(x) = (<f)(x)). It may well happen that at 
Tc there appear many different minima of T while at higher temperatures the only solution 
is M(x) = 0. If so, the presence of the source a would lift this degeneracy by selecting 
the closest solution M(x) ^ (i.e. having the largest overlap with a). The propagator 
Ghs{x — y) would then be the average of M(x)M(y) over all these solutions, that is over all 
possible fields a. A natural idea is thus to define a partition function Z solutions as the sum 
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of the Gibbs weig hts e~^ M ) of all the solutions of = 0. It would be very interesting 

to confirm or disprove that the equality between Z so i utions and the partition function of the 
field <p eqn.(||) also holds in the glassy phase below To- 
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